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We compare the failure probabilities of ensemble implementations of quantum algorithms which 
use pseudo-pure initial states, quantified by their polarization, to those of competing classical prob- 
abilistic algorithms. Specifically we consider a class algorithms which require only one bit to output 
the solution to problems. For large ensemble sizes, we present a general scheme to determine a criti- 
cal polarization beneath which the quantum algorithm fails with greater probability than its classical 
competitor. We apply this to the Deutsch-Jozsa algorithm and show that the critical polarization 
is 86.6%. 



I. INTRODUCTION 

There are two general paradigms for implementing 
quantum algorithms [1]. In the first, the quantum algo- 
rithm is implemented on a single quantum system with 
the appropriate number of qubits and which can be pre- 
pared in a suitable pure state and is amenable to projec- 
tive measurements. Most quantum algorithms are writ- 
ten with this in mind. In the second paradigm, the al- 
gorithm is implemented on an ensemble of of identical, 
non-interacting quantum computers. This is the situa- 
tion with conventional room temperature, solution state 
NMR implementations, in which case the ensemble con- 
sists of approximately 10^° molecules [2-8]. 

In ensemble implementations each ensemble member 
undergoes the same unitary evolution as its companions 
and algorithms for the two paradigms arc typically most 
similar in this respect. However, they differ in the initial- 
ization and measurement stages. In general an ensemble 
quantum computer can only be prepared in a mixed state, 
so that the state of any single ensemble member is not 
known with certainty. Also, the output from an ensemble 
quantum computer is an average of individual ensemble 
member measurement outcomes. The initialization and 
measurement issues have led to modifications of quantum 
algorithms for ensemble realizations. 

The conventional approach to ensemble quantum com- 
puting initializes the ensemble in a pseudo-pure state, 
for which various preparation techniques have been pro- 
posed [2, 3, 9-11] and which has the form 
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where n is the number of qubits, \tpi) is a known pure 
state and < e < 1 is called the polarization. The idea 
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is that under the collection of unitaries required to im- 
plement a quantum algorithm, L^aig, the density operator 
transforms to 
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and this is followed by measuring the expectation value 
of a traceless observable. The identity component of of 
Pfinai does not contribute to this measurement outcome 
and it is as though the pure state algorithm represented 
by Uidg\ipi) has been implemented. 

Much of the discussion of ensemble quantum comput- 
ing on pseudo-pure states has focused on the scaling 
properties of the polarization with respect to the prob- 
lem's input size [12] or the presence of entanglement in 
these [13]. In particular, most pseudo-pure state prepara- 
tion schemes result in polarizations which diminish expo- 
nentially as the number of qubits increases, thus result- 
ing in exponentially decreasing output signal strength. 
However, a promising new approach using NMR with 
parahydrogen induced polarization attains high polariza- 
tions and appears to avoid these problems [14]. 

Here we consider how well an ensemble quantum al- 
gorithm, for a given polarization and ensemble size, per- 
forms in relation to competing classical probabilistic algo- 
rithms. We propose a criterion, considering the ensemble 
size as one of the resources, for which an ensemble algo- 
rithm can be compared fairly to a classical competitor. 
We then use this to ask, for a certain class of problems, 
whether there is a critical polarization below which the 
quantum algorithm fails with greater probability than 
the classical algorithm. 

The remainder of this paper is organized as follows. 
In section II we provide a general scheme for compar- 
ing the performance of ensemble quantum algorithms to 
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their classical counterparts. We only consider algorithms 
for which the output is obtained after measuring a sin- 
gle qubit. In section III we apply the general scheme to 
the Deutsch-Jozsa algorithm determine the critical po- 
larization below which the quantum algorithm fails with 
greater probability than a classical random algorithm. 
Finally, the appendices contain much of the mathemati- 
cal derivations of various essential results. 



II. PERFORMANCE OF ENSEMBLE 
QUANTUM ALGORITHMS VS. CLASSICAL 
PROBABILISTIC ALGORITHMS 

We consider problems which take one of many possible 
inputs and determine into which of two possible classes 
the input falls. Any classical algorithm to solve one of 
these could be designed to write the output to one bit; 
those inputs returning "0" fall into "class 0" , and those 
returning "1" fall into "case 1." We assume that a quan- 
tum algorithm exists, which, when applied to a collection 
of qubits in an appropriate pure initial state, can deter- 
mine the input class with certainty. It is convenient to 
split the collection of qubits into a single qubit target reg- 
ister, on which a measurement will reveal the input type, 
and a remaining n-qubit argument and workspace regis- 
ter as may be required by the algorithm. This quantum 
analog proceeds as: 
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for "class 1" 
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where the subscripts denote the argument/ workspace 
and target registers and \4>o)a and |(/)i)a are normalized 
but not necessarily orthogonal argument register states. 
The input class is revealed following a computational ba- 
sis measurement on the target qubit. 

On an ensemble quantum computer initially in the 
pseudo-pure state of Eq. (1), the typical protocol [2, 15] 
for determining the input class is based on the expecta- 
tion value for the target qubit. 



for "class 0" 
for "class 1". 



(5) 



This evidently allows one to distinguish the input class 
by "measuring an expectation value" (provided that the 
polarization is suitably large for detection in a particu- 
lar experimental setup) and checking whether it is +e or 
— e. However, for an ensemble with a finite number of 
members AI and whose final state is mixed as in Eq. (4), 
the random nature of the target qubit measurement out- 
comes on individual ensemble members generates statisti- 
cal fluctuations which will yield outcomes that are almost 
never precisely {(7z)t = ±£- It is then essential to elabo- 
rate the protocol for deciding the input class, determine 
the probability with which this gives a correct result and 
compare this to a classical probabilistic algorithm which 
uses the same resources. 



The protocol which we advocate replaces {az)t by a 
suitable sample average of computational basis measure- 
ment outcomes over all the ensemble members. We as- 
sume that a computational basis measurement is per- 
formed on each ensemble member and that each mea- 
surement outcome is scaled to be compatible with the 
eigenvalues of cr^, i.e. let Zj = +1, Zj = —1 correspond 
to the outcome of the measurements associated with pro- 
jectors Pq — |0)(0| and Pi = respectively. These 
yield a sample average 
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which typically approximates {0^)^ well as Af — > 00. This 
leads to the decision protocol: 
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input is "class 0" , 

guess the input class with probability 
1/2 for either type, and 
input is "class 1" . 



(7) 



This amounts a majority vote on the number of individ- 
ual ensemble member outcomes which are Zj = -1-1 or 
= — 1 or a completely unbiased guess whenever the 
numbers of the two outcomes are identical. Let M+ be 
the number of times that that Zi = +1 and M_ the num- 
ber of times that Zi = —1. It is straightforward to verify 
that 



AM 



M 



where AM := A/+ — M_ represents the excess of positive 
measurement outcomes. The protocol of Eq. (7) assumes 
the best possible resolution in the measuring apparatus. 
That is, one can distinguish between AM = ±1 (for M 
odd) or AM = -2,0, or 2 (for M even). We refer to 
this as the best resolution case. We shall later general- 
ize this to arbitrary measurement resolution and demon- 
strate that the best resolution case is optimal. 

The probability with which the quantum algorithm 
misidentifies the input type can be determined by con- 
sidering the various routes to failure. The probabil- 
ity that that a "class 0" input will be misidentified as 
"class 1" will be denoted as Pfaiibesto a-rid the prob- 
ability that a "class 1" input will be misidentified as 
"class 0" as pfaii best 1 • Assuming that an input is cho- 
sen from "class 0" with the same probability as from 
"class 1," the quantum failure probability is Pfaiibost'' 
(Pfaii best o+Pfaii best i)/2. Now supposc that the algorithm 
is run with a "class 0" input. The input will be misiden- 
tified if M+ < M_ or if an incorrect class is guessed when 
— A/_ . The probabilities with which these occur can 
be derived from those for measurement outcomes on in- 
dividual ensemble members. In this case it follows from 
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Eqs. (2) and (4) that 
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Similarly if the algorithm is run with a "class 1" input 
the failure probability can be determined by switching 
M_i_ with M_ in the conditions for misidentification and 
Zi = +1 with Zi = —1 in Eq. (8). The symmetry in 
these situations implies that Pfaiibcsti = Pfaiibosto and 
thus Pfaiibcst*^ ~ Pfaiibcsto- Sincc measurements on each 
ensemble member amount to a Bernoulli trial the "class 
0" failure probability is a cumulative binomial distribu- 
tion. The precise form of this depends on whether M is 
even or odd. For odd M, the case M+ = M_ cannot 
occur and 



The best resolution case assumes that the measure- 
ment apparatus allows one to distinguish between two 
circumstances where the values of AM differ by as lit- 
tle as 2 and thus values of z which differ by as little as 
2 /M. In a general resolution case we assume that one can 
only distinguish between two situations where the values 
of AM differ by a resolution of at least R, which could 
depend on A/. In the context of the protocol of Eq. (7) 
this means that outcomes for which — i?/2 < AM < R/2 
can be regarded as pure noise. The maximum magni- 
tude of the sample average associated with this noise is 
\z\ = R/2M and noting that the maximum sample aver- 
age associated with any outcome has magnitude \z\ = 1, 
the signal to noise ratio is represented by R/2M. This 
can be used as a guide to precise behavior of the resolu- 
tion as a function of ensemble size, which may depend on 
the details of the apparatus. Regardless of these details, 
the decision protocol for the general resolution case is: 
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indicating the dependence of the failure probability on 
polarization and ensemble size. For even M, the case 
Mj^ = Af_ can occur and 



Pfaiibest'(£> M) = Pr(A/_ > A/+) + - Pr(A/_ = M+) 
= Pr(A/_ > AL + 1) + 1 Pi.(J\/_ = i^) 

= E 
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z > 7^ ^ input is "class 0" , 

guess the input class with proba- 
bility 1/2 for either type, and 
input is "class 1". 

(11) 

Note that the best resolution case is represented by 
R = 2. The symmetry in this protocol again results in 
Pfaii'^ ~ Pfailo • The "class 0" input failure probabilities are 
more conveniently expressed in terms of Af_. To do so, 
note that unequivocal failure, i.e. z < —5^1 corresponds 
to AM < -\R/2'] and, since 2A/_ = M - AM this is 
equivalent to Af_ > \{M + [i?/2])/2]. For convenience 
define the minimum number of occurrences of Zi = —1 
needed for unequivocal failure as 



Ar„ 



\{M+\R/2^)/2^. 



(12) 



Clearly Afmin > Af/2. Also, it is easily shown that the 
ambiguous outcome ^jj > z > — 2T7 equivalent to 
A/min - 1 > Ar_ > Af - Afinin + 1- Thus thc quantum 
algorithm fails with probability 



Pf^il^(£, Af, AAin) = Pr(A/_ > AAin) + - Pr(AAin - 1 > A/_ > M - M^un + 1) 



i (Pr(7\/_ > A/,„i„) + Pr(A/_ > A/ - A/,„i„ + 1) 
1 
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Several important properties of this general quantum fail- fixed M and Ai,jii,i, Pf^^-^{£i M, Afjnin) is a monotonically 
ure probability are proved in Appendix A. First, for 
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decreasing function of e and 



Pf,/(l,M,M„,i„) = 0. 



(14) 
(15) 



The former corresponds to a maximally mixed initial 
state, for which the algorithm produces a maximally 
mixed final state and any decisions about input classes 
amount to unbiased guesses. The latter case corresponds 
to a pure initial state, for which the algorithm never fails. 
Second, for fixed e and M, as the resolution decreases, 
i.e. M„iin increases, Pf^i^{e,M,M^i^) increases. Thus 
the best resolution case provides a lower bound on the 
failure probability for the quantum algorithm, as is to 
be expected. This bounding property is important since 
it appears to be easier to arrive at certain results for 
the best resolution case than the general resolution case. 
Two important results regarding the best resolution case 
are also proved in Appendix A. First, if M is odd then 
the best resolution case failure probabilities for M and 
M + 1 are equal. Second, if M is odd then the best res- 
olution failure probability for Af + 2 is strictly less than 
that for M unless e = or e = 1 (both statements re- 
quire fixed e). Thus, in the best resolution case at least, 
it is advantageous to using ensembles of increasing size. 

In general there are no closed form expressions for cu- 
mulative binomial distributions of the sort encountered 
in Eqs (9), (10) and (13). However, the following result 
due to Bahadur [16] can give good approximations, par- 
ticularly for M^oo. IfO<p<l,m and n are positive 
integers, and 



i?„(m) := J2 



Ap (I-P) 



(16) 



then, provided that np ^ m ^ n 
np{l — p) 



where 

An{m) : 



1 



(m — np)^ 



Bn{m) Anim) (17) 



{m + l){l~p) 
(m + 1) — (n + l)p 



(18) 

Consider first the best resolution case, in which case it 
is only necessary to consider situations where M is odd. 
It is straightforward to verify that the conditions for Ba- 
hadur's approximation are satisfied for the cumulative 
binomial distribution of Eq. (9). The factor on the left 
side of Eq. (17) becomes 



np{l-p) 
(m — np)'^ 



= 1 



(l-e2 



(1/VM + VMe) 



(19) 



and thus tends to 1 as M ^ oo provided that vAfe 
oo as M — > oo (this will be shown to applicable to the 
Dcutsch-Jozsa algorithm). In such cases the quantum 



error probability is well approximated by Eq. (18) after 
the correct substitutions for m, n and p. Now consider the 
general resolution case. Bahadur's approximation applies 
to the first term on the right of Eq. (13) since M^i^ > 
M/2 but in general the conditions are not satisfied for 
the second term on the right of Eq. (13). In other cases 
it is shown in Appendix A that it applies to the second 
term on the right of Eq. (13) when e > [fl/A/. Thus 
provided that R scales as RqM" where Rq is constant 
and < a < 1, the approximation applies for almost all 
e as Af ^ oo. The result analogous to that of Eq. (19) 
must be determined for each term on the right of Eq. (13). 
For M 1 the first term gives 



1 + 



np{l-p) 



1 



(l-e^) 



{m-npf ■ {\R~\/2^/M + ^/MeY 
while for the second term it gives 



(20) 
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np{l~p) 
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(l/Vl/ - \R~\ /2 VA7 + ^/Mef 

(21) 

Again, these tend to 1 as as Af —* oo provided that 
VAfe — *■ cxD as Af — *■ oo and the quantum failure prob- 
ability is well approximated using Eq. (18) twice with 
appropriate to, n and p. 

It remains to compare the failure probability for a 
quantum algorithm to that for competing classical proba- 
bilistic algorithms. This is easiest for algorithms, such as 
the Deutsch-Jozsa algorithm or search algorithms, which 
solve problems with the aid of an oracle. In these the 
input is a function / drawn from one of two classes. The 
only aid allowed is an oracle which can evaluate / at 
any possible argument. The task is to determine the in- 
put type with the fewest oracle queries. We henceforth 
restrict the discussion to such oracle query algorithms. 
We are concerned with cases where M is very large since 
these are typical in NMR realizations and also the quan- 
tum failure probability in the best resolution case de- 
creases as M increases. However, the ensemble size must 
be included in the count of resources and we do so by 
incorporating this into the total number of oracle queries 
(this has been used in the context of ensemble realizations 
of the Deutsch-Jozsa algorithm on thermal equilibrium- 
type states[17]). Suppose that Uaig invokes the oracle q 
times. Since f/aig is applied to each ensemble member, 
the aggregate number of oracle queries is Q := Mq. Thus 
a quantum algorithm using q queries per quantum com- 
puter operating on an ensemble with A/ members must 
be compared to a classical probabilistic algorithm which 
uses Q oracle queries. Denote the classical failure prob- 
ability with Q oracle queries by Pf^i\(,Q)- It is assumed 
that < Pf^{\{Q) < 1/2 and that Pfs^i\{Q) decreases as 
Q increases. Then the critical polarization is the min- 
imum £ required for the quantum failure probability to 
drop beneath the classical failure probability, is obtained 
by solving Pf,ii=(0) = Pi^^'{Mq) = p^^^{e, M, M^i^) for 
e. Since Pf^j['^(e, Af, Afmin) decreases monotonically from 
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1/2 to with increasing e, there will be a unique critical 
polarization, e(M), for each M. 

The precise behavior of e(Af ) depends on the behav- 
ior of the ratio the quantum failure probability to the 
classical failure probability as a function of M as well as 
the behavior of the resolution as a function of AI. This 
is somewhat simplified by considering the best resolution 
case since it bounds the quantum failure probability for 
the general resolution case from below and will provide a 
lower bound on e(M). Thus consider the best resolution 
case. If the critical polarization is bounded from below 
in the sense that there exists Mq and eo > such that for 
M > Mq, e(M) > Eg then the conditions for Bahadur's 
approximation apply and it gives (see Appendix B) 



The circuit for the standard Dcutsch-Jozsa quantum 
algorithm is illustrated in Fig. 1 where the gate opera- 
tions are defined on computational basis states as 

^i-^) = 7fE(-ir'|y) (25) 

y=o 

for the Hadamard gate and 

Uf\x)\y)^\x)\y(Bf{x)) (26) 

for the oracle. These are extended linearly to arbitrary 
superpositions of quantum states. It is straightforward 



e{M) = ^l-[M{p,^^HMq))^Y/'' 



(22) 



for large M. 

For example, consider a classical probabilistic algo- 
rithm for which Pfaii'^(Q) ~ l/c'^ where c > 1. It is shown 
in Appendix B that if 7\/ > 2/ logc then e > ^1 - l/c^. 
This satisfies the conditions leading to Eq. (22) and gives 
a critical polarization in the best resolution case of 



In the asymptotic limit, M^/^' ^ 1 as A/ ^ oo and 



(23) 



e(M) 



1 



(24) 



The general resolution case depends on the behavior 
of the resolution i? as a function of M. However, if R 
scales as RqM" where Rq is constant and < a < 1, 
then Bahadur's approximation again applies and it is 
straightforward to show that as M — > oo, A/min ^ Af/2 
which approaches the best resolution case. It follows that 
Eqs. (22) to (24) apply to this situation as well. 



III. EXAMPLE: THE DEUTSCH-JOZSA 
ALGORITHM 

The Deutsch-Jozsa problem [18] considers functions 
/ : {0, 1}" {0, 1} which are guaranteed to be either 
constant or balanced. A balanced function yields for 
precisely half of the = 2" possible arguments and 1 
for the remaining half. The task is to identify the func- 
tion type using the minimum number of invocations of an 
oracle which can evaluate f{x) at any x ~ 0, . . . N — 1. 
The approaches for determining the function type with 
certainty are well-known [18, 19]; classically, in the worst 
case, the function must be evaluated for 2"~^ -I- 1 differ- 
ent arguments; if two different inputs have yield different 
outputs it is balanced but if all inputs return the same 
output it is constant. 
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FIG. 1: Quantum circuit for the standard version of the 
Deutsch-Jozsa algorithm. The actions of the gates are defined 
in the text and the algorithm terminates with a computational 
basis measurement on all n control register qubits. 



to demonstrate that if / is constant then final state of 
the two registers is [V'finai) = |0...0)|1). while, if / is 
balanced, [^Afinai) = Yl,x=i Notably, for a bal- 

anced function, the state |0 . . . 0) does not appear in the 
argument register superposition. Thus an ri-qubit com- 
putational basis measurement on the argument register 
reveals the function type. This quantum algorithm re- 
quires just one oracle invocation to accomplish this (giv- 
ing q = 1). 

In the language developed earlier, the constant func- 
tions correspond to "class 0" and balanced functions to 
"class 1" and the algorithm should be modified so as to 
yield a single bit output. This is accomplished by an ad- 
ditional multiply-controlled NOT as illustrated in Fig. 2. 
If / is constant, the final state is of both registers is 
IV'finai) = |0o)|O)7 while if / is balanced, the final state 
will be IV'finai) — for some (irrelevant) \(t>j). Thus 

a computational basis measurement on the target qubit 
reveals the function type. Note that the extra multiply- 
controlled NOT gate can be decomposed into a sequence 
of 0{n^) basic one and two qubit gates [20]. 

The framework developed earlier can be used to com- 
pare the performance of ensemble realizations of this al- 
gorithm to its classical probabilistic counterparts. The 
classical probabilistic algorithm proceeds by evaluating 
/ on M < N/2 + 1 distinct arguments. If all outputs are 
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register 
(n qubits) 



Target 
register 



|0) 



|0) 



|1) 




FIG. 2: Modified quantum circuit which produces a single bit 
output for the Deutsch-Jozsa problem. The final gate is a 
multiply-controlled NOT which applies a NOT to the target 
register when every argument register qubit is in state |0). 



the same / is identified as constant, whereas if two out- 
puts differ / will be identified as balanced. This can only 
fail when a balanced function happens to return the same 
output for all M arguments. Assuming that a balanced 
or constant function is chosen with equal probability, it 
is shown in Appendix C that the probability with this 
occurs is well approximated by 



1 

2M 



(27) 



provided that M < N/2. 

The critical polarization is determined by solving 

Pf.^{^^M) = ^. (28) 

For the best resolution case, the approximation of 
Eq. (22) with c = 2 gives 



£(A/) = Wl 



1 



(A/) 



l/M 



(29) 



We note that a better approximation for intermediate 
ensemble sizes is 



(30) 



These are illustrated, along with data obtained by numer- 
ically solving Eq. (28), in Fig. 3. In the limit M oo, 
Eq. (24) imphes e yfifl = 0.866025. By comparison a 
standard room-temperature, solution state NMR realiza- 
tion on 500 MHz spectrometer, using pulsed pseudo-pure 
preparation schemes [2, 21, 22] typically has e < 10~^. 
A more promising but more complicated method [14] us- 
ing parahydrogen induced polarization to produce a two 
qubit ensemble quantum computer has attained s = 0.9. 
It should be noted that, to date, all NMR realizations of 
the Deutsch-Jozsa algorithm [4, 23-29] have had n < 5 
and M ~ 10^° ^ N/2 and, by our criteria, a classical al- 
gorithm with comparable resources would determine the 
function type with certainty and thus outperform these 
realizations. 
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0.8' 
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FIG. 3: Critical polarization vs ensemble size for the Deutsch- 
Jozsa algorithm. The solid line is generated via Eq. (30) while 
the dashed line is generated via Eq. (29). The squares display 
data obtained by solving Eq. (28) numerically for the best 
resolution case, while the asterisks are for a resolution R = 



IV. CONCLUSION 

In conclusion, we have provided a method for com- 
paring the performance of ensemble versions of quantum 
algorithms whose output is extracted from a measure- 
ment on a single qubit to their classical probabilistic 
counterparts. We have applied this to realizations of the 
Deutsch-Jozsa algorithm and calculated the minimum 
polarization required for the quantum algorithm to out- 
perform the classical probabilistic algorithm. Our calcu- 
lations indicate that the standard room temperature so- 
lution state NMR approach attains polarizations several 
orders of magnitude too small but that newer approaches 
using parahydrogen induced polarization attain suitable 
polarizations for the ensemble quantum computer to out- 
perform the classical probabilistic algorithm. 



APPENDIX A: QUANTUM FAILURE 
PROBABILITY 



The following useful representation of cumulative bi- 
nomial distributions [30] can be verified by repeated in- 
tegration by parts: 



k—m 

where Ip{x, y) is the incomplete beta function defined by: 
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1. Behavior with respect to e 

It is trivial to show by direct substitution that 
Pfaii'^ll' Mmin) = 0. Now consider 



Pf,,''(0,Af,AW) = i J2 

fc=A/„i„ 
M 

k=M-M„,i„ + l 
1 



A/\ / 1 \ / 1 \ 



2MH 



fc=M„ 



and, since Mmin > M/2, 



1 ^ 
7+T X] 



fe = A/-A/n,in + l 



1 

'^(0,Af,M„,i„)= 2M E 



fe=M„i„ 



A/ 
A: 

Jl/mi„-l 



fc=M-M„i„ + l 

It is straightforward to show that 



E 

fc=A/-J\f„i„ + l 



2 E 

fc=L(A/+l)/2j 



and thus 

Pf^il'»(0,A/,A/„,in) 



1 

2^ 
1 

2' 



A/ 



E 

A;=[(A/+1)/2J 



(A3) 



Now consider the behavior as e increases. We show 
that 



de 



< 



(A4) 



for < e < 1. To prove this, note that derivatives of cu- 
mulative binomial distributions are easily computed us- 
ing Eqs. (Al) and (A2); 



dBn{m) d 
dp dp 



Ip{m, n — m + 1) 



T{n + 1) 



r(m)r(?i — m + 1) 



p"'-\i-py 



> 



provided that < p < 1. Equation (13) shows that 
the quantum failure probability is just the sum of two 
positively weighted cumulative binomial distributions 
with p = "'^i ^-nd applying the chain rule proves the 
result. 



2. Behavior with respect to Afmin 



We show that, for any fixed e and M, 



Pf^i,'^(e, A/, Af,„i„ + 1) > Pf^/(£, M, A//^in) (A5) 



provided that M/2 < M^i^ < A/ - 1. Let 



^PM Pfaii''(£> M, A/„in + 1) - Pf,ii'^(e, M, M„ 



Then 
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Since < e < 1, the term between brackets is positive. since Ip{x +1,2:)+ Ip{x, x + 1) = 2Ip{x, x). This proves 
This proves the result. the result. 



3. Best resolution case behavior with respect to A/ 

Consider the quantum failure probability for the best 
resolution case when M is odd. We show that the quan- 
tum failure probability remains constant if one additional 
ensemble member is added: 



ftailbe.t'(e: ^ + 1) = Pf,,bc.t'(£, M) . (A6) 



Let 



^PM Pfaiibcst'(e. + 1) - Pfaiibcst'(e, M). 
Then Eqs. (9), (10) and (Al) imply, with p := (1 - e)/2, 



M+l M+1 



+ ll (M±l + 1 M+l) 



T (M+l M+l" 

^pV 2 ' 2 , 



Now consider passing from M to M + 2. We show that, 
for the best resolution case and M odd. 



Pfailbe.t'^(£: M + 2) < ft,, M), (A8) 



with equality only when e = or e = 1. Let 



^PM := ftailbest'(e' M + 2) - ft,iibest'(£> M). 



(A7) -jj^gj^ Eqg^ g^jj^j (Aj^) jj^^piy^ ^^^-^ p--(i- e)/2, 
I 



^Pm{p) 



A/+1 



M 

M+l 

— 2~/ "'O 



1 ALhi^ 
2 ; 



^(M-i)/2(^ - t)(^^-i)/2 [2t(i _ t)(M + 2) 
I 



i^] dt. 



Now consider 



, M \ M+l M+l 

9{P) ■■=-[ M+l (l-2p)p 2 (1-p) 2 . 



APPENDIX B: QUANTUM FAILURE 
PROBABILITY VS CLASSICAL FAILURE 
PROBABILITY 

1. Applying Bahadur's approximation 



It is straightforward to show that 
dp dp 

and that ^(O) = Apm(O). Since both functions are con- 
tinuous it follows that they are identical. 



/ M \ M+l M+l 

APM = -(^jU±lj(l-2p)p 2 (l-p) 2 , 



However 0<p<l/2for0<£<l, and thus ApM < 0. 
For e = and £ = 1, corresponding to p = 1/2 and p ^ 
respectively, A.pM — 0. This proves the result. 



Consider a circumstance where it is known that there 
exist a positive integer Mq and £o > such that for 
M > Mq, e{M) > Eq and where the resolution scales as 
RqM" where i?o > is constant and < a < 1. Then 
\fMe oo as A/ — * cx). Thus the term 



1 



npjl - p) 
(m — np)'^ 



in Eqs. (19), (20) and (21) tends to 1 as Af ^ oo. It 
remains to approximate An{m) in Eq. (18). Clearly for 
the resolution which scales as described above, Afmin ~ 
A//2 if A/ > 1 and Eq. (13) imphes 



1 



1 



ft,ii'^(£, A//, Af„i„) ~ - AMiM/2) + - AmUM + l)/2) 
~Am{M/2). 
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For M :$> 1 this gives 



M/2 



l + e 



{m/2) ( 2 
/ M \ jl-s^y'/^ (l+e) 
\M/2) 2^1 2e ■ 



M/2 



(Af/2 + l)(l + £)/2 
(A//2 + 1) - (A/ + l)(l-e)/2 



(Bl) 



The binomial coefficient can be approximated using Stir- 
Ung's formula 



2. Exponential classical failure probability 



for n > 1. Thus 

M 
M/2 



Ml 



{M/2)\{M/2)\ 



^/2^ M*^ 



ttM {Ml2Y'i 



Consider the case ViB:i\{Q) = l/c*^ where c > 1 and 
Q = M(7 is the total number of oracle queries over the 
entire ensemble. Wc shall prove that there exist a positive 
integer Mq and eo > such that for M > Mq, e(M) > eq- 
The strategy is to consider the ratio of the best resolution 
case quantum failure probability to the classical failure 
probability, 



givmg 



The critical polarization is determined 
Pfa/(e'^^'^^min) ^Pi^^{'{Mq) and this yields 



1 



p,,^%Mq)V2^I 



1 



2/M 



by 



(B2) 



Now consider M ^ 1. Since we have assumed that 
£(M) > eo 7^ 0, the factors on the right which are con- 
stants or contain e are approximately 1. Thus 



p,,,CiMq)VMj , 



giving Eq. (22). 

Note that this can be improved for intermediate sized 
M by retaining the factor of \/27r in Eq. (B2). However, 
we found an even better approximation for the case of 
the Deutsch-Jozsa algorithm by solving numerically for 
e(M) for small M and substituting in to the last fraction 
on the first line of Eq. (Bl). This explains Eq. (30). 



Pfail ratio (£,A^) - 



Pfailbcst''(£.^) 



and to show that for some Mq and eo > 0, 
Pfaii ratio (e,M) > 1 wheu M > Mq and eo > e. 
This establishes that the critical polarization, for which 
Pfaii ratio (e, M) = 1 , is bouudcd from below by eo and 
this applies regardless of the resolution, since the best 
resolution case provides a lower bound for polarization. 

The crux is to establish that for sufficiently small e, 
Pfaii ratio (e, -M) iucrcascs as M increases. Note that for 
odd M, pfaii ratio (e,M + 1) > Pfaii ratio (e,M) for any e 
since the best resolution case quantum failure probability 
remains constant while the classical failure probability 
decreases. Thus consider 



Ap(e) := Pfailratio(£, M + 2) - Pfail ratio (e, M) 



for odd M. Then using Eq. (Al) with p := (1 - e)/2, 



J 



A/(e,M)=c«^^ [c2%( 



M 

M+l 



A/+1 



M+l M+l 



)] 



{t{l - t))(^^-i)/2 ^2t(t - 1){M + 2y - ^] 



dt. 
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Then 



dAf 
de 



M+l 
2 



2M 



r 



which is negative if 



< Jl _ M+l 29 
^ ^ ' M+2 



(B3) 



Since the right hand side of Eq. (B3) increases as M 
increases, < when e < -y/l — 2c^^i/3. 

Now consider Pfaiiratio(0, M) = c'^^' /2. Then 
Pfaiiratio(0,Af) > 1 whcn M > Mq := [log2/qlogc] > 
1. But there is some polarization e' > such 
that Pfaiiratio(eVMo) = 1. Then choosing Eq = 
minje', y/1 — 2c^'^i/3} ^ (note that this is indepen- 
dent of M) imphes that Pfaii ratio (e, M) > 1 for M > Mq 
and e < Eq. Finally this implies that e{M) > Eq ^ for 

M > Mq. 



anccd function yields M identical outputs. If the choice 
of balanced functions is unbiased, then the probability 
with which this occurs is the number of balanced func- 
tions for which the first M arguments all return the same 
result divided by the total number of balanced functions. 
The number of balanced functions which return (or 
equivalently 1) for the first M arguments is 



and the total number of balanced functions is 



N-M 
\N/2-Mj 
M \ 



KM/2/ 



Thus the probability of misidentifying a balanced func- 
tion is 2 y (7V/2) where the factor of 2 counts 
both the cases which output and those that output 1. 
For ^ 2M, this can be approximated using Stirling's 
formula. Thus the classical failure probability is 



APPENDIX C: CLASSICAL FAILURE 
PROBABILITY FOR THE DEUTSCH-JOZSA 
PROBLEM 

In the classical probabilistic algorithm, / is evaluated 
on M < N/2 arguments. The algorithm fails when a bal- 



N - A-I ' 
N/2 - M 



N 

N/2 



{N - M)\[N/2)\ 
N\{N/2-M)\ 



'{N-M)N/2 fN/2-M 
N{N/2 - M) \ N-M 



M<JV/2 1 

' 2^' 



M 



iN-M)iN/2) \ 
{N/2-M)N J 



N 



(CI) 



Thus the classical failure probability tends to pbai2/2^^ 
where pbai is the probability with which a balanced func- 
tion (vs a constant function) is chosen. In the unbiased 
case considered in this paper, pbai = 1/2. Thus the clas- 
sical failure probability is well approximated by 1/2^^ 
provided that M < N/2. 
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